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Abstract 

The analytical approximations for the moderating neutrons flux density like Fermi 
spectra, widely used in reactor physics, involve the probability function for moderating 
neutron to avoid the resonant absorption obtained using some restrictive assumptions 
regarding the acceptable resonances width. By means of multiplicative integral (Volterra 
integral) theory for a commutative algebra an analytical expression for the probability 
function is obtained rigorously without any restrictive assumptions. 

1 Introduction 

The compound- nuclear reactions may be divided into resonant and non resonant. It is 
known that the nucleus energy may possess only a discrete set of values corresponding to 
its energy levels. However the idea of the levels with strictly fixed energy is true only for 
the ground levels (non-excited levels) of the stable nuclei. All the other (excited) levels 
do not possess the fixed energy - they are more or less smeared over energy. A width Tj 
(i is the level index) of such smearing can be estimated using the uncertainty relation for 
time and energy. According to this estimate, AEi = Tj = h/ri, where Tj is the lifetime for 
current level. The nucleus may be excited through the interaction with some projectile 
particle only if this excitation energy induced by interaction corresponds to the energy 
levels spacing. Therefore a compound nucleus may form only in the case that the projectile 
particle energy is such that the excitation energy of the compound nucleus induced by 
this interaction fits into the corresponding level's interval of uncertainty Tj (if the nucleus 
had been in its ground state prior to interaction). If the width of compound nucleus 
energy levels is less than the distance between them and the projectile particles possess 
some fixed energy, the reaction may take place only trough some i th single level. The 
reaction cross-section dependence on projectile particle energy would display a resonant 
behavior. Correspondingly, this kind of reactions is called resonant. If the levels are 
located so densely that the distances between them is less than their width, they merge 
together and the reaction may take place for any projectile particle energy. This kind of 
reactions is called non-resonant. The form of resonances is known to be described by the 
Breit-Wigner function. 

As is known, the resonance region of the reactor neutron nuclear reactions of fissile 
nuclides is located, depending on neutrons energy, in the 0.5-^1000 eV range (e.g. [3] and 
Fig. 1). 

Determining the neutrons spectrum for a moderating medium with resonant absorption 
remains an important task for numerous problems of nuclear reactor physics. Numerical 
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Figure 1: 235 JJ capture and fission cross-sections dependence on neutrons energy. Dashed line 
represents the neutrons scattering cross-section dependence. 

methods of integro-differential equation solution for the moderating neutrons flux density 
in the medium with resonant absorption are used as well as the widespread analytical ap- 
proximations for the flux density like Fermi spectra. Fermi spectra involve the probability 
function for moderating neutron to avoid the resonant absorption. The specific expres- 
sion for this probability function is usually obtained using some restrictive assumptions 
regarding the acceptable resonances width. The narrow, broad (infinite mass approxima- 
tion) and intermediate resonance approximations are used for that purpose. Therefore 
depending on the properties of nuclei composing the moderating medium with absorption, 
one of the above approximations is used. 

We show in the present paper that the probability function for a moderating neutron 
to avoid the resonant absorption (let us note that the neutron absorption is the opposite 
event, and therefore the absorption probability function is equal to one minus the proba- 
bility function to avoid the absorption) may be written down in a form of multiplicative 
integral (Volterra integral) [1]. The analytical expression for this function was rigorously 
obtained using the multiplicative integral theory for a commutative algebra [2] without 
any restricting assumptions. 

2 Probability function for avoiding the resonant 
capture with narrow resonance approximation (stan- 
dard approach) 

Fermi spectrum for a moderating medium with resonant absorption, according to [3, 4], 
has the form: 

$*(E) = =£-.<p{E), (1) 

where ip(E) is the probability function for a moderating neutron to avoid the resonant 
absorption; S is the total volume neutron generation rate; C = S (&^s) /^s, £i is the 

i 

average logarithmic decrement of energy loss; Eg is the macroscopic scattering cross- 
section and is the macroscopic absorption cross-section of the i th nuclide; E| = ^ + 

i 

S* is the total macroscopic cross-section of the fissile material, = ^5 i s the total 

i 

macroscopic scattering cross-section of the fissile material. 
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Figure 2: Two distant resonances E\ and E2 {AEx, AE2 are the widths of the first and second 
resonances respectively). 



A following standard approach is used in reactor physics [3, 4] in order to obtain the 
probability function <p(E). According to the standard slowing-down theory and Fig. 2, 
which represents the resonance region for neutron absorption cross-section in a simpli- 
fied visual manner, the probability for neutron to avoid the absorption within the first 
resonance region may be written in the form: 



1 



AEi 



where S a is the macroscopic absorption cross-section; E\ is the first resonance energy and 
AE\ is its width (see Fig. 2). 

The probability for neutron to avoid the absorption within the second resonance region 
with energy Ei and width AE2 has the following form: 



¥>2 



1 



AEi 



£(£ s + £ a ) Ex 



1 



AEo 



£(x 5 + £ a ) E 2 



(3) 



Thereby by analogy we have the probability function for a moderating neutron to 
avoid the resonant absorption: 
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(4) 
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where A is a number of resonances. 

In order to pass from multiplication to summation, we take a logarithm of the function 
<p(E) defined in (4): 
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Inip(E) = JJln I 1 



where x 
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Next every multiplicand in expression (5) is expanded into series in terms of a small 

AEi 



quantity x% with an assumption that 
obtain the following expression: 

In </>(£) w 



<C 1. Limiting to the first expansion term, we 
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(6) 



Let us divide the resonance region into m energy intervals of width AEj and assume 
T, a = in the intervals between the resonances. Thus for (6) we obtain: 
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Passing from the logarithm to an exponent and making m approach the infinity, for 
entire resonance region from Eq to Et we have: 



ip(E) ~ lim exp I — N 



exp 



\ So 



(8) 



The Fermi spectrum (1) with the probability function (p(E) in exponential form (8) is 
called a narrow resonance approximation. 



3 Multiplicative integral and its application to a 
commutative algebra A 

Basic structures related to a multiplicative integral occur in different areas of mathemat- 
ics, mechanics, physics and - as it follows from the sections 1 and 2 above - in the nuclear 
reactions and nuclear reactor physics as well. However one cannot say that the multi- 
plicative integral theory is currently "particularly popular" among mathematicians and 
physicists. Let us define multiplicative integral according to [1]. Let A be an arbitrary 
associative topological algebra with an identity E, and f(t) is a function of real variable 
t which takes on the values in the algebra A, a, b G R, a ^ b, [a, b] = {t G R, a ^ t ^ b}, 
T is the partition of an interval [a, b] with the points to = a, ti, ^2> • • ■ , £n-l> *n = b, 
ti < t i+ i,i = 0, 1,2, ... ,n - 1; l(T) = max(t i+ i - t{), At { = t i+ i - U. 
Let us consider the following product: 

II (/. T) = {E + /(to) Ato) x (E + /(ti) Ah) x...x(E + f(t n )At n ) . (9) 

If for any variation T such that l(T) — > the product Y\ (/> T) converges to some limit, 
this limit is called a multiplicative integral of the function /(t) over the interval [a, b] and 
is denoted by 

u b 

jE + f(t)dt. (10) 
By analogy, for f(t) and T one can derive the following: 



n II (/' T) = (E + f(t n )At n ) x(E + f(t n ^)At n ^) x...x(E + /(t )At ) . (11) 

If for any variation T such that Z(T) — > the product n n (/, T) converges to some 
limit, this limit gives another type of the multiplicative integral for the function f(t) over 
the interval [a, b] and is denoted by 



j E + f(t)dt. (12) 



The multiplicative integrals (10) and (11) are called direct and backward respectively. 
Next let us consider a particular case of multiplicative integral application. Let the 
algebra A be commutative. Then, according to [1, 2], 
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E + f(t)dt = lim (E + f(t )At )x(E + f(t 1 )At 1 )x...x(E + f(t n )At n ) 

J Z(T)->0 



lim 

Z(T)->0 



£ + Yl f( f i) At i + f^fit^AtiAtj + ^ f{ti)f{ti)f(t k )AtiAtjAt k + ... 

i<j i<j<k 

(13) 



= lim [cr + <7i + ci2 + . . . + a n ] , 

where o"o, o\, 02, . . . , <x n are the elementary symmetric polynomials in variables f(to)Ato, 
/(ti)Ati, f(h)At 2 , .., f(t n )At n . Let 

PO) Pi; P2; • • •) Pn be symmetric Newton's polyno- 
mials in the same variables: 



Pk = Yl [/(*i) A *il* » fe = l s 2,. 



,n. 



(14) 



It is easy to see that 



lim = < 

i(T)-X) 



0. 
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fc > 0,A;^ 1 



(15) 



where j f(t)dt is the Riemann integral. 

a 

Expanding the polynomials cr^ in terms of pi (k,l = 1,2, ... ,n) and applying the 
equality (15), we calculate 



lim Ufc. 
l(T)-*o 



(16) 



It is known that 



p m - p m -\0\ + Pm-2<?2 + ■ ■ ■ + Pl<T m „i (- l) m 1 + (-l)"V m = 0, TO = 1, 2, . . . , n, (17) 

whence it follows that 



limcri = limpi = J f(t)dt, 
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f(t)dt 



lim OV 



n! 



f(t)dt 



(18) 



and as a result (taking into account that s x = ^ + % + ^X + ^\ J r 
n — > oo when l(T) — > 0) we obtain 



•• + ^r + 



. and that 



J E + = exp (| /(t)di j . 



(19) 



Due to Eq.(13) symmetry, for the backward multiplicative integral we, obviously, ob- 
tain the following equality: 



n b 

J E + f(t)dt = exp ( / IH)dl J . (20) 

a 

4 Multiplicative integral and its application to 
the derivation of a probability function for a mod- 
erating neutron to avoid the resonant absorption 

According to the stated above, it is easy to establish the correspondence between the 
theory of a probability function for a moderating neutron to avoid the resonant absorption 
which is a real scalar function (and therefore a particular case of commutative algebra 
A) and a multiplicative integral. Indeed, according to the section 3, the identity E of 

this commutative algebra A is the real unit 1, and f(E) = — = is a scalar 

£(£5 + £ a )£ 

function of a real argument E which takes on the values in the algebra A, Eo,Ef £ R, 
E < Ef, [E ,Ef] = {E eR,E ^ E < E f }, T is the partition of an interval [E Q ,E f ] 
with the points to = E ,h,t2, ■ ■ ■ ,t n -i,t n = Ef, U < U+i,i = 0, 1,2, ... ,n - 1; l(T) = 
m&x(t i+1 - ti), Aij = ti+i - U. 

Then for the probability function to avoid the resonant absorption, according to (4) 
with the partition T of an interval [Eo,Et] we obtain: 

nE f 

- ( s \ n r 

tp(E) = lim TT 1 - = AEi = lim TT (1 + /(E 4 )A£<) = / E+f(t)dt, 

Eo 

(21) 

where f(E) = —= . 

£(£ s + £ a )£ 

Consequently, taking into account (20), we derive the final expression: 



(p{E) = exp 



Sa dE 



. £(Es + £ a )£ 
\ E J 



(22) 



5 Conclusion 

In the present paper we showed that the probability function for a moderating neutron to 
avoid resonant absorption may be written down in a form of a multiplicative integral. Us- 
ing the results of the multiplicative integral theory application for a commutative algebra 
the analytical expression for the probability function for a moderating neutron to avoid 
the resonant absorption was rigorously obtained. The expression for Fermi spectrum in- 
corporating the obtained probability function turned out to coincide with the expression 
for Fermi spectrum known as the narrow resonance approximation. However, the former 
one is an exact expression and therefore does not have any restrictions. 

We should also note here that since the narrow resonance approximation is used for 
another types of nuclear reactions, the obtained results may be generalized for a broader 
range of nuclear physics problems. 
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